VISCOSITY SOLUTIONS AND AMERICAN OPTION PRICING 
IN A STOCHASTIC VOLATILITY MODEL OF THE 
ORNSTEIN-UHLENBECK TYPE 



ALEXANDRE F. ROCH 

Abstract. In this paper, we study the valuation of American type deriva- 
tives in the stochastic volatility model of Barndorff-Nielsen and Shephard [2]. 
We characterize the value of such derivatives as the unique viscosity solution 
of an integral-partial differential equation when the payoff function satisfies 
a Lipschitz condition. 



1. Introduction 

In their seminal paper, Barndorff-Nielsen and Shephard |4J introduced a model 
that has been shown to describe particularly well financial assets for which log- 
returns have heavy tail distributions and display long range dependence. In this 
model, the volatility of the asset is described by an Ornstein-Uhlenbeck type pro- 
cess with a pure jump Levy process acting as the background driving process. An 
empirical study was made in [1] and showed from exchange rate data that suitable 
distributions for the Levy process are the so-called generalized inverse gaussian dis- 
tributions from which well understood examples are the normal inverse gaussian 
(studied in and the gamma distribution. 

The BNS model has been studied from different points of view. Benth et al. 
[5] considered the problem of optimal portfolio selection. Nicolato and Vernados 
[TO] have studied European option pricing and described the set of equivalent 
martingale measures under this model. To evaluate these types of options, the 
authors propose the transform-based method and a simple Monte Carlo method. 

In this paper, we consider the pricing of American options with the use of 
integral-partial differential equations (IPDE). Although our technique can be sim- 
plified and used for European options and certain path dependent options such 
as barrier options (see [6j for a definition and examples), we will mainly concen- 
trate on American type derivatives which have not been studied for this model. 
The main difficulty in this case is the lack of Lipschitz continuity of some of the 
coefficients of the IPDE. 

The connection between viscosity solutions of IPDE's and Levy processes has 
been studied in the literature by various authors. Pham [11) considered a general 
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stopping time problem of a controlled jump diffusion processes. However, his re- 
sults do not apply here because the Lipschitz condition on the coefficients is not 
satisfied in our current setting. Cont and Voltchkova [6] studied barrier options 
and Barles et al. [2] established the connection between viscosity solutions and 
backward stochastic differential equations. In these papers, the stock price consid- 
ered is modelled by a stochastic differential equation with jumps driven by a Levy 
process. The main difference between the BNS model and these models is the 
presence of stochastic volatility. However, we will see that the lack of smoothness 
of the solution to our IPDE will also lead us to consider the notion of viscosity 
solutions as presented in [7]. 

The rest of the paper is organized as follows. In Section[2l we present the model 
and recall the results of Nicolato and Vernados [10] regarding the set of equivalent 
martingale measures. Section [3] is devoted to the continuity of the value function. 
In Section 2J we prove that the value function is the viscosity solution of the 
associated IPDE and the uniqueness of the solution is presented in Section [5l 

2. Levy Processes and the BNS Model 

Let T > 0. We consider the stochastic volatility model of Barndorff-Nielsen and 
Shephard [4] for the price process of an asset, denoted by S = {S t }o<t<T and 
defined on a filtered probability space (fl^F, {^Ft}o<t<T, P)- We thus assume that 
the log-return X t = log(S t ) of the asset satisfies the following stochastic differential 
equation: 

(2.1) dX t = {p + f3V t )dt + ^V t dB t +pdZ xt 
with 

(2.2) dV t = -XVtdt + dZxt 

in which /^,/3 S K, A > and p < 0. B = {B t }Q<t<T is a Brownian motion and 
Z = {Z t }o<t<T is the background driving Levy process (BDLP) under the physical 
measure P. In this model, Z has no gaussian component and the increments are 
positive. Z and B are assumed to be independent, and F = {!F t }o<t<T is the usual 
filtration generated by the pair (B,Z). The positivity of the jumps of Z insure 
that the process V is always positive. We denote by W the Levy measure of Z. 

Suppose Q is a probabilty measure equivalent to P under which S is a martin- 
gale. We are interested in American-type derivatives of the form 

U t = esssup TeT ^ T y t ^[e' r ^-^h{X T )\T t ] 

in which h is the payoff function and Tt is the set of all stopping times with values 
less or equal to T. Since {X t }o<t<T and {V t }o<t<T are Markov processes, U t can 
be written as a function of (x, v, t), say 

U t = u(x,v,t)= sup E Q (e- rT h(X*' v )) 

in which (Xf' v ) t >o is the process X for which X$ — x and Vq — v. We also denote 
by (V t ") t >o the process V starting at Vq = v at t = 0. 
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2.1. Equivalent Martingale Measures. We start by summarizing the results of 
Nicolato and Vernados [10] concerning the set of equivalent martingale measures. 
In order to do so, we define the set 

y = ly : [0,oo) -> [0,oo) ; J (Vvi x ) ~ lfw(x)dx < oo 

and M.' as the set of all equivalent martingale measures Q such that Z is still 
a Levy process under Q independent of J5, possibly with a different marginal 
distribution. 

As in jlOj . we impose the following conditions on the process Z: 

(CI) The process Z is given by the characteristic triplet (0,0,1^) so that the 
cumulant transform is given by 



/>OC 

k(9) = log{E[exp(0Zi)]} = / (e 

Jo 



l)W(dz); 



for values of 6 for which this expression is defined. 
(C2) 6 = sup{0 G R | k(6) < oo} > 0; 
(C3) lim k(6») = oo. 

Remark 2.1. Assumption (C2) implies that there exists 9q > such that 

(e ?oZ - l)VF(dz) < oo. 



For 2 > andn > I, we have < z n < 2±{e e ° z -l), so that /i„ := J°° z n M^((iz) < 

oo. Furthermore, Assumption (C2) is a sufficient condition for the process Z to 
have finite moments of all orders. 

The following theorem was proved in [lOj . 
Theorem 2.2. For all Qe M , there exists y € y such that 

dX t = (r-XK y (p)-^V^\dt+^V t dBf + pdZ Xt , 

/>oo 

in which k v (8) = / (e 8x — l)y(x)w(x)dx, 



and = B t ~ J (y/Vs)^ 1 ^ — fj, — (j3 + ^)V S — \n y (p))ds and Z\ t are respectively 
a Brownian motion and a Levy process under Q. w y (x) — y(x)w(x) is the Levy 
density of Z\ under Q and k v {9) is the cumulant function. 

In the remaining part of this paper, all expectations will be with respect to a 
chosen EMM Q, unless specified otherwise, and W and B will denote the associated 
Levy measure and the Brownian motion associated to Q. 

Let = lx R + x [0, T) and assume for a moment that u is Lipschitz in (x, v) 
and 

(2.3) ueC 2 ' hl {0), 
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that is u is differentiable with respect to v and t, and twice differentiable with 
respect to x. We can then apply Ito's formula to U to find 

du du — 

(2.4) dU t = (— + C[u])dt + —VVtdBt + dV u 

in which 

I du du I d^u 

C[u] = (r--v-Xny(p) + Xp^)--X(v-^)- + -v^ 

f°° du du 

+A J (u(x + pz,v + z, t) — u(xv, t) — (p z -q^ + z ~Q^))W{dz). 

and Vt is the Q-martingale given by 

dV t = J (u(X t -+pz,Vt-+z,t)-u(X t -,Vt-,t) 

- (pz?^(Xt- , Vt- , t) + z^(X t - , Vt- , t))) N{dz, Xdt) 

du du \ — 

pz— {X t _ , V t - , t) + z—(X t _ , y t _ , t) j N{dz, Xdt) 

in which N(dz, dt) = N(dz, dt) — W(dz)dt and N(dz, dt) is the random measure of 
the process Z. Since J ^y/VtdB t is a Q-martingale, if it can be shown that e~ rt Ut 
is also a martingale we can then expect u to satisfy the following integral-partial 
differential equation (IPDE) 

du 

— (x, u, i) + £[w](ar, v, t) — ru(x, v,t) = 

if u(x, > Otherwise w(x,v,t) = and this IPDE can be written as 

du 

(2.5) max( — (x, v,t) + C[u](x, v, t) — ruCx,v,t),h(x) — u(x, v,t)) = 0. 

dt 

It is clear also that the function satisfies 

(2.6) u[x, v, t) = h(x) for v = or t = T. 

Condition 12.31 is in fact very restrictive and most of the time not satisfied. 



Despite this problem, we will see that u can still be regarded as a solution of this 
equation in a weaker sense. 

3. Continuity of the Value Function 
Recall the definition of the value of an American option with payoff h: 

(3.1) u{x,v,t)= sup E (e- rr h{X^ v )) . 

In the rest of this paper, we will assume that h is positive and satisfies the Lipschitz 
condition, in other words 3K > such that V(xi, X2) S M 2 

(3.2) \h{x x ) - h{x 2 )\ < K\xi - x 2 \. 

For instance, the payoff function for an American put with strike X > is h(x) = 
max(JT — exp(x), 0) and satisfies this condition. 
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Our goal is to show that the function u satisfies the IPDE (|2.5p in some weak 
sense. In order to give meaning to this IPDE for a function u that doesn't sat- 
isfy basic differentiability conditions, we introduce the idea of viscosity solutions 
following Crandall and Lions [8]. Let W be the set of functions / : O — > R that 
satisfy 

\f(x,v,t) -f(x',v',t)\ w . . 

sup t \ , < oo Vt e 0,T . 

(x,t,),(x',t ; ')eRxR + l + \x-x'\ + \v-v'\ 

Definition 3.1. The function u £ C (O) DW is a viscosity subsolution (superso- 
lution) of |JDJ)-|]LIP ifV(x, v,t)eO and V^eWn C 2 ^\0) such that 

(i) ip(x,v,t) = u(x, v, t) and 

(ii) V(ar', v', f) e O tp(x', v', t') > u(x', v', t') (<), 
then 

max ^^(avM) + £^ ( Xj ^, t) - rV>(a;, M) ; 

(3.3) %) - u(a;, u, i)) > (<), 

(3.4) and u(x, v, t) = h(x) for v — or t = T. 

The function u is a viscosity solution if it is both subsolution and super solution. 

Remark 3.2. As noted in [S] (p. 317) the condition ^ S W is sufficient to have a 
well defined integral term in C[tp]. In fact if ip € Wfl C ' ' i/ien 

/ °° (^(x + pz, v + z, t) - i>{x, v, t) - (pz%(x, v, t) + z%(x, v, *))) W(dz) 



< / Cz'W(dz)+ \ C(l + \z\)W(dz) <oo 

'2<»7 



/or any 77 > 0. 



An important property of viscosity solutions is the continuity of the function. 
It is the content of the following proposition. 



Proposition 3.3. When h satisfies the Lipschitz condition the function u 

is continuous and in W. 

Proof. In this proof, we will assume for simplicity that r = 0. The generalization 
to r > is straightforward. Throughout, C is a positive constant that can change 
from line to line. 

We start by showing the continuity of u with respect to (x,v), uniformly in t. 
We have the following representation of the volatility process: 

Y, r = »e~ xt + [ e- Xs dZ 



t = vt. -t- / e azj\ s . 
Jo 

We define the integrated variance process started with Vq = v by 

vr = fv-ds. 

Jo 
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By Equation (12. 2j) . we find that 



V?dt 



-<-dV t v + dZ 
A 



\i ) 



so that we have the following representation of the integrated variance process 



(3.5) 
(3.6) 

in which e(t) 



V 



A 



A 



dZ 



Xa 



(I 



ve{t)+ / e(s)dZ X s- 

10 



l-e" 



We also have the following identities: 





— x — 




AV t 


■= v? 


- V t v = Ave~ xt , 


av; 


■■= vf 


* - V t v '* = Ave(t), 


AX t 


■= Xt 


— X t = x — x — 




:= Ax 


~ l -Ave(t)+Mr\ 



vy' - Jv? dB s , 



with Ax — 



and Av = v ' — v. 



Using the Lipschitz condition on h, we obtain 

\u(x', v', t) — u(x, v, t)\ — 



sup Eh(X* < v ) - sup Eh{X*> v ) 

T£T T _ t ' T E T'Y — t 



< sup E 

reT T -t 



< C sup E 

rer T _ t 



h{X% ' v ) - h{Xp v ) 



Then, 



\u(x', v', t) — u(x, V, t)\ < C (\Ax\ + \Av\+ sup E\M?' V | 



T£T T _ t 

Letting Q = ct({Z s }o< s <t), the cr-field generated by the BDLP Z up to the 
maturity T, we find that {M t ' }t>o is a Q V J^-martingale. Thus, {\M^' V \}t>o is 
a Q V J^-submartingale and Doob's theorem applies. In other words, 

sup E\M?< V '\ < E( sup e(\M?' v '\ 



< E 



( E (|M^;i|a))< Je(e((m^; 
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Also, 

E ((M^;) 2 |e) 



= £ \v s v ' -2y/v a "'V s v + V 8 v 



)ds 



T-t 



Ave'^ds + 2 



T-t 



v v 



\J (V s v ) 2 + V s v Ave- Xs ds 



- r- 

Jo 



Ave~ Xs ds + 2 



T-t 



-V?Ave- Xs 



V s v + sj {Vg) 2 + VgAve~ Xs 



ds 



< 



/ 3|Au| 
Jo 



e- Xs ds = 2,\Av\e{T - t) < Z\Av\T. 



And thus we proved the continuity of u in (x, v) uniformly in t since 
\u(x',v',t)-u(x,v,t)\ < C (\Ax\ + \Av\ + ^/lA^j . 
In particular u £ W because of the following inequality 

\u(x',v', t) — u(x, v, t)\ < C (\Ax\ + \Av\ + y/\A$) 
< 2C(l + |Aa;| + |Av|). 
The next step of the proof is to show 

E sup \XZ> V -X? V \^Q and E sup \V S V - V t v \ -» 

t<s<t' t<s<t' 

as |i — t'| — *• 0. This is easily obtained by first observing that 

Esu Pt<s<t , \x^-xr\ 



< -E sup \V S V >*-V t v '*\+-E sup 

* t<s<t' t<s<t' 



VjdB y 



+ pE sup \Z Xs 

t<s<t' 



< - w | e (t')-e(t)|+E 



e -Kt's) dz 



As 



-CJE \V t V* - V t v '*\ + P E \Z xt , - Z xt \ 



< - <t)\ + (1 + P)B \Z\v - Zxt\ 



-v\e(t')-e(t)\+E\Z Xt ,-Z xt \. 
As for the process V, 

E sup \V S V -V t v \ < |1 — e -A(*'-t) | | -h E 



t<s<t' 



f e-^'-°Uz Xs 



< \l-e- x «'-V\E\V t v \+-E\Z Xt ,-Z xt \, 
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Since V t v < v + Z XT for all t < T, 

E sup \V S V - V?\ < C(v + EZ XT )\t' - t\ + E \Z xt , -Z Xt \, 

t<s<t' 

and we need to show that E \Z\t' — Z\ t \ — > when \t' — t\ — ► 0. 

We mentioned earlier that condition (C2) implies that the moments of Z t are 
finite for all orders. Thus Z is uniformly integrable. Since Z is also continuous in 
probability, it is continuous in C\ and the conclusion follows. 

Let's now show continuity with respect to time. Let < t < t' < T. Take 
r e T T - t and define r' = r A (T - t'). Then, 

E,(e~ rT h(X^ v )) 

= E L- rT ' h{X*l v )\ + E (e- rT h{X*' v ) - e- rT ' h(X*; v )^j 
< u(x, v, t') + E (e- rT h(X*' v ) - e- rT 'h(X*; v fj . 
From this inequality, we readily find that 

\u(x,v,t')-u(x,v,t)\ < E sup \Xf> v - XZ'O.t,] 

T-t'<s<T-t 

which converges to zero as \t — t'\ — * 0. 

Global continuity follows from the following inequality 

\u(x',v',t)—u(x,v,t)\ < \u(x , v , t') — u(x, v, t')\ + \u(x,v,t')—u(x,v,t)\ 

and the fact that the first bound is independent oft'. □ 

4. Viscosity Solutions 

This section is devoted to the viscosity solution property of the value function 
u. In order to prove that u is a viscosity solution of (12. 5|) . we need the following 
dynamic programming principle. It is a consequence of the martingale property 
of the Snell envelope stopped before its optimal stopping time and it is the key 
property needed in the proof of the subsolution property. 

Lemma 4.1. Let e > 0, (x,v,t) € O and define the stopping time 

r e = inf{0 < s < T - 1 1 e- rs u(X*< v , V s v ,t + s) - e < e~ rs h(X*- v )}. 
Then, 

(4.1) u(x,v,t) =E[e- rT 'u(X*l v ,V^,t + T e )]. 

Proof. For some constant C, we have that 

nE(\l WXr) > n} h(X T )\) <E(h(X T ) 2 ) <C + CE(X 2 T ) 

for all r <E T T . We know that X T = X + tt + V* + \/V^dB s + pZ Xr and that 
< V* < Y£ < \(V + Z XT ) from Equation As a result, X 2 T < 4(X + 

rT) 2 +4±(V + Z XT ) 2 +4(f ( ;VV s dB s ) 2 +4p 2 Zl T < C + CZ 2 T + C (J^ VV s dB s ) 2 
for some constant C large enough. Hence ~EX 2 < C + CEZ 2 T + CEVf < oo for 
all r S Tt- As a consequence, sup Tg7 - T E (|l{ /l (x T )> rl }/i(A r )|) converges to as n 
grows to infinity, i.e. the collection {e~ rT h(X T ) : r G Tt} is uniformly integrable. 
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Hence we find that the process (e~ rs h(X s ))o< s <T is of Class D and we can apply 
the results of [9] to get the result. □ 
The proof of the solution property of u makes use of the following lemma. 

Lemma 4.2. Let t < T and e > 0. Suppose u(x,v,t) — h(x) > e. Then Q(r e < 
s) — ► when s — > 0. 



Proof. Let i] > e such that r\ < u(x, v, t) — h(x). 

First we show that e~ rT 'u{X*i v , V£ , r e ) - e - rT ' h(X?l v ) < e almost surely. 
For some sequence s n [ r e , e~ rSn u(X^ v , Vg n ,s n ) < e~ rSrl h(X^ v ) + e for n large 
enough. In this case, since (Xf^ 3 , Vg n ) converges to (X*f, V^i) in C 1 , we can take 
a subsequence if necessary and find that \u(X%' v , , s n ) — u{X*f , V^i, r £ )| — > 
and h(Xf' v ) — > h{X*f) a.s. with n — > oo. Taking the limit, we find 

e-^ M (X^,y; e ,T £ ) = lim e- r8 »uffi»,^ )Sn ) 

n — udo 

*h(X*i v ) + e a.s. 



< lim e 

n — >oc 



Since u is continuous with respect to t, we find that r\ < e rs u(x,v,t + s) 
e~ rs h(x) for s small enough. Then, for s small enough, 
Q(r e < s) 



< 
< 
< 



i(e- rTC (u(x,v,r^ - h(x)) + e~^\ h {X x T ?) - u{X x T f ,V^)) > V - e) 

H 



U{X,V,T*)-U{X X T ?,VZ,T*) 



h{x) 



V?. - v 



>S 2 ) +q( 



X% 1 



>s. 



for some constants 5 2 > and S3 > 0. By the continuity in probability of the 
processes X et V, we know that this expression goes to zero when s — > 0. □ 
We can now show that u is a viscosity solution. 

Theorem 4.3. When h satisfies the Lipschitz condition H3.2\) , u is a viscosity 
solution oflPDE JO) . 



Proof. We already know that it is continuous and in W. 

Let's start by showing that u is a supersolution of ()2.5|) . Let (x, «, t) € and ip 
satisfy the conditions given in the above definition of supersolutions. By definition, 
for any At > 0, 



e 



-rAU 



E(u(X^,Vl t ,t + At))- U (x,v,t) 



> E (e-^(X^, Vl t ,t+ At) - V {X x ' v , \%,t)) 



E 



At 



+ J ■£(X^ 1 V:\t + s)e- rs ^dB s + 
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in which ^ x,v is the martingale defined by 



(VpCf + pz, +z,t + s)- ^{X^, V s v _ ,t + s) 
-z (M + ^\ (X^,V?_,t + s))N(dz,\ds) 



Since 



is also a martingale, we have the following inequality 

> jf * E (e- rs (-r^ + ^+ CW){X? », V?,t + *)) da, 

in other words, dividing by At and taking the limit as At — > 

> -rip(x,v,t) + ^-(x,v,t) +£[ij)](x,v,t). 
ot 

Since, by definition, u(x, v,t) > h(x), u satisfies Equation (|2.5p . To prove that u 
is a viscosity subsolution of (|2.5|) . let (x,v,t) £ O and ?/> satisfy the conditions of 
the above definition for subsolutions. If u(x,v,t) — h{x), the inequality (|3.3|) is 
satisfied. Otherwise, let 

< e < u, t) — ft,(a;). 
We know from Lemma 14.11 that 

= E( e -( At ^) U (^ ATE ,^ ATE ,t+(AtAr e )))- U (x, W ,^ 
< E(e- r ^'^(X^,Vl tAT „t+(AtAr^)) -#r, V ,t) 

(4.2) = E(j^\-^ { ^ + ^ +C [^X^,V s v ,t + S )d S \ 



for any At > 0. Knowing that Q(r e < At) — * when At — * by Lemma I4T21 
dividing the preceding inequality by At and taking the limit to 0, we get the 
desired result by Lebesgue's dominated convergence theorem. □ 



5. Comparison Principles and Uniqueness of the Solution 

In this section, we prove a comparison result from which we obtain the unique- 
ness of the solution of the IPDE. In proving comparison results for viscosity solu- 
tions, the notion of parabolic superjet and subjet as defined in Crandall et al. [8 
is particularly useful. Setting y = (x,v), we define the parabolic superjet and its 
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closure by 

J 2 ' + u(y,t) = |(p,g,i)eRxl 2 x5 2 such that u(y' \t') - u(y,t) < 



p(t'-t) + q- (y' -y) + ^(y 1 — y) T ■ A ■ (y' - y) 



+o(\t'-t\ + \y'-y\ 2 ) as (f, y') - (t, y)} 

2 + f 

J' u(y,t) = Hp, g, A)= lim (p„,g„,yl„) such that 

L n — >oo 

(p„,q„,A n ) e J 2 ' + u(y n ,t n ) and (y„, i„) -> (y, t) j. 

The subjet and its closure are then defined similarly by 

J 2 ^u(y,t) = -J 2 >+(-u)(y,t) and 

J 2 -u{y,t) = -J^ + (-u)(y,t). 

We then have the following lemma which is essentially proved in [2\ (lemma 
3.3). 

Lemma 5.1. // the function u S C°(K x M + x [0,T]) is a viscosity subsolution 
(resp. super solution) of \2. 5}) then \/(x,v,t) G R x R + x [0, T) and V(p, q, A) G 

2 + 2 — 

J ' u(x,v,t) (resp. J ' u(x,v,t)) 

(5.1)max(j> + C q ^ A [u, tp](x, v, t) — ru(x, v, t) ; h(x) — u(x, v, t)) > (<), 
in which 

Cf A [u^](x,V,t) = (r-^v-\ K y(p) + \p^)q^-\(v-^)q^ + ±vA u 

f £ dtb dib 

+A / (ip(x + pz,v + z, t) - ip(x, v, t) - (pz-- + z--))W(dz) 
J Q ox Ov 

f°° dip dtp 

+A / (u(x + pz, v + z, t) — u(x, v, t) — (pz— — 1- z—-))W(dz) 

ox ov 

for some ip £ C ' ' and < £ < 1. 

Pham [IT] obtains the uniqueness of the solution when the coefficients of C 
satisfy Lipschitz conditions on K 2 x [0, T]. For 8 > 0, define O 5 = fx (8, oo) x [0, T). 
Then, the coefficients of £ satisfy the Lipshitz conditions on O s and using the ideas 
of uniqueness proofs in the literature we can show a comparison principle on O s . 
This result will then be used to show the uniqueness of the solution on O. 

Theorem 5.2. Let e > 0, and u\ be a subsolution and u 2 a supersolution of 
on O s such that 

u\(x,v,t) < u 2 (x,v,t) +e 
for t = T or v = 8. Then u\(x, v, t) < u 2 (x, v, t) + ee r ( T ~'' for all (x, v, t) S O s . 
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Proof. AnIPDE of the form a ^ x Q f^ + C[ip](x,d,t)-rip(x,ti,t) = for (x, <!?,*) £ 
O and ip(x, i?, T) = h(x) was shown to have a unique solution in [2] when the 
coefficients of C satisfy some given Lipschitz conditions. In fact when (x, i?, t) and 
(I'.fi'le^we have | V¥ - V$\ < ^ \&' - i?| and so the operator £ satisfies the 
assumptions made in [2]. The extension of the uniqueness result to our current 
setting is straightforward and we only give the main details. 

We first show that u\ — u 2 is a subsolution of a related IPDE. Suppose ip £ 
W H C 2 and u\ — u 2 — ^ attains a maximum at (j/q, to) G C 5 - Set 

w(yi,U2,t,s) = wi(yi,t) - u 2 (y2,s) 

and 

(f>(yi,y 2 ,t,s) = ^\yi ~ y 2 \ 2 + - s\ 2 + tpiy!^). 

Since wi and u 2 are in W, the function w — <fr attains its maximum (y*, y 2 ,t* , s*) 
(which depends on e, a) in O s x O 5 . By a classical argument in the theory of 
viscosity solutions we can show that \\y{ — y 2 \ 2 , ^\t* — s*\ 2 — > when e, a — > 
and 

(2/*>2/2>**= s *) (yo,yo,to,t ) 

when e, a — > 0. 

Applying Theorem 8.3 of Crandall et al. [7] to the functions w and <f>, we find 
matrices Y\ , F2 such that 

(a+^iylXlb + D^ylX)^) £ J 2 ' + Ul {yt,t*) 

(-o,-6,-y 2 ) e J 2 ' + (- W i)(2/ 2 *, s *). 

with a = — a*) and 6 = -(y\ — y 2 ) and for < £ < 1 the inequalities 

mK („ +4 , [ „ I ,, 1(9 ,, ) + g«,, )+ „ : (s: ,, ) _„ lW ,, ); 

and 

mayi^a + C^ Y2 [u 2l -ip](y 2l s*) -ru 2 (y 2 ,s*) ; fcfca) ~ "2(2/21 s *)) <° 

are satisfied. Write these two expressions as max(A, B) > and max(C, D) < 0. 
Thenmax(A-C,B-£)) > 0. Now, B-D = h{x* x )-u x {yt, t*)-h{x* 2 )+u 2 {y^ a*), 
and because /i is Lipschitz |/i(x*) — /i(x2)| — > when e, a — > 0. Thus £> — D — > 
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U2(yo,to) ~ u>i(yo,to). On the other hand it was shown in [2] that 
A-C < r (u 2 (vl8*)-u l W,1*)) + ~<~-\ + ±M-vS\ 2 

+A J (il>(xl + pz, v{ + z, t*) - ip{yl,t*) - z(p, 1) ■ £>v) W(dz) 

+^ J (<t>{ x * +pz,v* + z, xl,v%,t*, s*) - 4>(x*,v* ,x* 2 , vl, t*,s*) 
-z(p,l)-(b + DiP(ylt*)))w(dz) 

J ($( x h v l, x 2 + P Z , V 2 +Z,t*,S*) - (j)(xl,vl,xl,V2,t*,S*) 

-z(p,l)-b 2 )w(dz), 



in which r' = (r — ^n y (p) + A/i 2 ). Using the fact that <b 6 W H C 2 we find letting 
£ — * and then e, a — > that 

(5.2) A — C < - r (ui{yo,h)-MvoM) + -£+ty- 

at 

Consequently, 

dib 

max(-r(ui - u 2 )(y ,t ) + -57(2/0, *o) + Cib(y ,t ), 

(5.3) -(«i-«2)(»d,*o))>0. 

As shown in [2] (see lemma 3.8), there exists a function x > 1 such that 

^ + £x-rx<0 

and for which the maximum 

M = sup ((ux - U3)(y, t) - /3 X (y, t))e r <'- T > 

RxB + x[ti,T] 

is attained at some point (yo^o)- Then 

(ui-u 2 -/3x)(|/,t) < ( M i- U2 -/3x)(2/o,to)e r(t °- t) . 
Let ?/>(y, i) = (3x(y, t) — {u\ — u 2 — Px)(yo, to)e r< - t °~ t \ Then ^> satisfies the properties 



in the subsolution definition, hence it satisfies Equation 15.31 But 

dib d\ 

(-^ + £ip)(y ,t ) = (/3-g+r(u 1 -U2-l3x) + l3Cx)(yo,to) 
< r(u 1 -u 2 )[yo,to). 

Hence, either (u± — U2)(yo,to) < 0, or vq — S or to = T and, in this case, (u± 
U2)(yo,to) < e by assumption. Hence, we conclude that 

( Ul -u 2 )(y,t) < /3 X (y,t)-/3 X (y ,i )e'-( t «- t ) + ( Ul - W2 )(2/o,to)e r(to - t) 
< p X (y,t) + ee r (*°-*>. 
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Sending /3 to zero we get ui < U2 + ee rlyT l > on R x {8, oo) x \t\, T}. As done in [2], 
we can repeat this argument as many times as needed to get u\ < U2 + ee r ( T ~*' 
on O 6 . □ 

A solution of (|2.5jl - (|2.6p is said to be minimal if it is less or equal to any other 
solution of (|2~5|) - (j2"l)|) . 



Theorem 5.3. u is the minimal viscosity solution of S2.5\) -s2 
Proof. Let S > and define 

u s (x,v,t) = sup E (e- rT h{X^ v )) 

TeT T -t-T<T 6 

in which 

t s = mi{s >0:V S V < 5}. 
Then u s is a viscosity solution of (|2.5p on O 5 with boundary conditions 

(5.4) u 5 (x,v,t) = h(x) for t = T or v = 6. 

The proof of this statement is essentially the same as the proof for the viscosity 
solution property of u. The main difference is that the maturity T is replaced by 
ts- Note that Vf > S for 5' > 5e XT , hence u s (x,v,t) = u(x,v,t) for all x € M, 
t < T and v > 5e XT . 

Let u be another viscosity solution of (|2.5p - (l2.6p . Then u is a viscosity solution 
of (12. 5p on O s with boundary values u(x, v, t) for t = T or v = S. Also, u(x, v, t) > 
h(x) = u (x, v, t) for t = T or v = S. By Theorem 15.21 we find that u > u s on 
O s . In particular, £t(:r,u, t) > u(x,v,t) for a; € K, t < T and u > <5e AT . Since J is 
arbitrary, u > u on C □ 

Following Pham [11], we denote by UC XtV (0) the set of functions defined on O 
uniformly continuous in (x, v), uniformly in t. We have already shown that the 
function u satisfies 

\u(x', v', t) - u(x, v,t)\<C (\x' -x\ + \v' -v\ + y/\v' - v\j . 

Hence, u £ UC x>v (0). Using the two previous theorems, we can show the unique- 
ness in UC x _ v (0). 

Theorem 5.4. u is the unique viscosity solution of I12.5\) - ![27B\) in UC x ,v{0). 

Proof. Let u E UC x , v (0) be another viscosity solution of (|2.5p - J2.6|) . Let e > 0. 
Then there exists S > such that < u(x, v, t) — u(x, 0, t) = u(x, v, t) — h(x) < e 
and < u(x,v,t) — u(x,0,t) — u(x,v,t) — h(x) < e for v < S. In particular, 
\u(x, v, t) — u{x, v,t)\ < e for all x, all t and v < S. Furthermore, by Theorem 15. 3[ 
we obtain that u(x, (5, t) < u(x, (5, t) < u(x, (5, t) + e, and u{x, v, T) = u(x, v, T) by 
definition. By the comparison principle of Theorem 15.21 we find that u(x,v,t) < 
u(x, v, t) < u(x, v, t) + te r{ - T ^ t ' ) for all (x, v, t) £ O s . Hence, u(x, v, t) < u(x, v, t) < 
u(x,v,t) + ee rT for all (x,v,t) £ O. Since e is arbitrary, we obtain the desired 
result. □ 
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